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Two-sublattice ferrimagnet, with spin-si operators Sn at the sublattice A site and spin-S2 op- 
erators S2i at the sublattice B site, is considered. The magnon of the system, the transversal 
fluctuation of the total magnetization, is a complicate mixture of the transversal fluctuations of the 
sublattice A and B spins. As a result, the magnons' fluctuations suppress in a different way the 
magnetic orders of the A and B sublattices and one obtains two phases. At low temperature (0, T*) 
the magnetic orders of the A and B spins contribute to the magnetization of the system, while at the 
high temperature {T* ,Tm), the magnetic order of the spins with a weaker intra-sublattice exchange 
is suppressed by magnon fluctuations, and only the spins with stronger intra-sublattice exchange has 
non-zero spontaneous magnetization. The T* transition is a transition between two spin-ordered 
phases in contrast to the transition from spin-ordered state to disordered state (Tiv-transition) . 
There is no additional symmetry breaking, and the Goldstone boson has a ferromagnetic dispersion 
in both phases. A modified spin-wave theory is developed to describe the two phases. All known 
Neel's anomalous M{T) curves are reproduced, in particular that with "compensation point". The 
theoretical curves are compared with experimental ones for sulpho-spinel MnCr2S4,-xSex and rare 
earth iron garnets. 

PACS numbers: 75.50.Bb, 75.30.Ds, 75.60.Ej, 75.50.-y 



I. INTRODUCTION 

The notions of ferrimagnetism and ferrimagnetic ma- 
terials were introduced by L. Neel [l[ for materials which 
spontaneous magnetization is a resultant of two or more 
components of non-parallel magnetic moments. Using 
a molecular field theory L. Neel predicted the nature 
of magnetization M versus temperature T. Anomalous 
M{T) curves arise due to the fact that each of the mag- 
netic moments approaches its own saturation value as 
a different function of temperature. In general there is 
nothing to limit the number of components. The simplest 
model consists of two alternating sublattices of unequal 
and antiparallel moments, with three molecular field co- 
efficients utilized to describe the exchange field effects: 
one ferromagnetic coefficient for each sublattice and a 
third for antiferromagnetic interaction. In many cases, 
the experimental curves are used to determine the molec- 
ular field coefficients. The most striking feature of the 
anomalous M{T) curves is the possibility "of compen- 
sation point" , a temperature Tc at which the magnetic 
moments of two sublattices are equal and opposite, so 
that the magnetization of the system is equal to zero 
M(Tc) = 0. The phenomenon of ferrimagnetism has been 
subject of extensive experimental investigation since its 
discovery. The phenomenological Neel's standpoint has 
been confirmed by many authors (see the review articles 

[a, 4, 5, 6:? ]). 

The earliest application of spin waves to ferrimagnet 
was made by H. Kaplan Q to calculate the dispersion 
LUk for the spinel-type ferrite. For small wave vector k he 
found a quadratic relation ujk = -Dk^, where D is a con- 
stant, as for ferromagnetism. To simplify the calculations 
Kaplan neglected the intra-sublattice exchange compared 
with inter-sublattice one, so that the calculations do not 
correspond to any real ferrimagnet. 



In the present paper I consider two-sublattice ferrimag- 
net, with spin-si operators Sii at the sublattice A site 
and spin-S2 operators 821 at the sublattice B site. The 
Hamiltonian of the system is 

H = — Ji Sii ■ Sij — J2 ■ 

+ JY,^u-S2, (I.l) 

where the sums are over all sites of a three-dimensional 
cubic lattice: denotes the sum over the nearest 

neighbors, <C «, j 3>a denotes the sum over the sites of A 
sublattice, <C i,j denotes the sum over the sites of 
B sublattice. The first two terms describe the ferromag- 
netic Heisenberg intra-sublattice exchange Ji > 0, J2 > 
0, while the third term describes the inter-sublattice ex- 
change which is antiferromagnetic J > 0. 

The true magnons of a two-spin system arc transver- 
sal fluctuations of the total magnetization which includes 
both the magnetization of sublattice A spins and magne- 
tization of sublattice B spins. The magnon excitation is a 
complicate mixture of the transversal fluctuations of the 
sublattice A and B spins. As a result the magnons' fluc- 
tuations suppress in a different way the magnetic orders 
on the different sublattices and one obtains two phases. 
At low temperature (0,r*) the magnetic orders of the 
A and B spins contribute to the magnetization of the 
system, while at the high temperature (T*, T/v) the mag- 
netic order of the spins with a weaker intra-sublattice ex- 
change is suppressed by magnon fluctuations and only the 
spins with stronger intra-sublattice exchange has non- 
zero spontaneous magnetization. At first sight the result 
that there is a temperature interval where a magnetic 
order is formed only by spins on one of the sublattices 
seems to be counterintuitive because the moment on the 
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one of the sublatticcs build an effective magnetic field, 
which due to inter-sublattice exchange interaction leads 
to finite magnetization on the other sublattice. This is 
true in the classical limit. In the quantum case the spin 
wave fluctuations suppress the A and B magnetic orders 
at different temperatures T* and Tn as a result of differ- 
ent interaction of magnons with sublattices' spins. The 
T* transition is a transition between two spin-ordered 
phases in contrast to the transition from spin-ordered 
state to disordered state (T/v-transition). There is no 
additional symmetry breaking and the Goldstone boson 
has a ferromagnetic dispersion in both phases. 

A modified spin-wave theory is developed to describe 
the two phases. By means of this method of calculation 
the thermal variation of magnetization is calculated for 
three cases: (si > S2) Ji » J2, (si > S2) J2 ^ Ji and 
(si = S2) Ji ^ J2- AH known Neel's anomalous M{T) 
curves are reproduced including the curve with a famous 
"compensation point" (second case). 

An important issue is the experimental detection of 
T* transition. In the case when T* is the temperature 
at which itinerant electrons start to form magnetic order 
the transition demonstrates itself through the change in 
the T dependence of resistivity. It is well known that the 
onset of magnetism in itinerant systems is accompanied 
with strong anomaly in resistivity Q. One expects the 
same phenomena when the itinerant electrons in ferri- 
magnets form magnetic order. There are experimental 
results which support the above interpretation of the T* 
transition. 

The paper is organized as follows: In Sec. II, a spin 
wave theory of the model Eq. p.l[) is presented. It is 
shown the non-adequacy of the customary spin- wave the- 
ory for description of the high temperature phase. Sec. 
Ill is devoted to the development of a modified spin wave 
theory. The calculations are accomplished for different 
choices of the model's parameters. The theoretical M(T) 
curves are compared with experimental ones for sulpho- 
spinel MnCr2S4-xSex and rare earth iron garnets Q 
and satisfying coincidence is obtained. A summary in 
Sec. IV concludes the paper. 



when the sites j are from sublattice A, and 



S+ = 51,. + ^Sl = -b+ J2s2 bp, 



S2j = ^2^ - ^Sl, = - y/2s2 - b+b, b, (II.2) 
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-S2 + b+b, 



when the sites j are from sublattice B. The operators 
Oj , aj and bj satisfy the Bose commutation rela- 
tions. In terms of the Bose operators and keeping only 
the quadratic terms the effective Hamiltonian Eq. ()I.l|) 
adopts the form 

H = siJi ^ {afai + apj - api - a.f'^j) 

+ S2J2 J2 i^th + bp, - bp, - bp,) (II.3) 

+ ^ [siJbpj + S2Japi - J^siS2 {ap'j' + aibj)] 
iv) 

To proceed one rewrites the Hamiltonian in momentum 
space representation 

4ak + ei bpk - Ik {apk + 6+afe) ] , 

feG-B,. 

(11.4) 

where the wave vector k runs over the reduced first Bril- 
louin zone B^ of a cubic lattice. The dispersions are given 
by equalities 

el = 4si Ji Ek + 6s2 J 

(II.5) 



el = 4s2 J2 Efc + 6si J 



with 



and 



efe = 6 — cos{kx + ky) — cos(fc2, — ky) — cos{kx + k^) 
cos{kx — kz) — cos{ky + kz) — cos(fc.y — k^) (II.6) 

7fc = 2 J -^si S2 (cos kx -I- cos ky + cos k^) (II.7) 

To diagonalize the Hamiltonian one introduces new 
Bose fields a^, a^, (3k^ (3^ by means of the transforma- 
tion 



II. SPIN- WAVE THEORY 

To study a theory with Hamiltonian Eq. (jl.ip it is con- 
venient to introduce Holstein-Primakoff representation 
for the spin operators 



2.S1 



2s- 



ap, (H.l) 



si — ap, 



Qk = Ukak + Vk Pt flfe = "fc "fe + '"k Pk 

(11.8) 

bk ^ UkPk + Vk bl ^ Ukf3^ + Vk ak 

where the coefficients of the transformation Uk and Vk 
are real function of the wave vector k 



Uk 



\ 




{et + ei)^~4^l 



Vk = signi'jk) 



\ 




(11.9) 
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The transformed Hamiltonian adopts the form 

H = Y, (e^ c.+a, + eP P+(3u + E^) , (11.10) 



with new dispersfons 



E? = I 



^{ei + sir ' 47I -si + st 



(11.11) 



el 



and vacuum energy 



El^l 



- El. - 



(11.12) 



For aU values of k 



H > 14 - ell 



and the dispersions are nonnegative E^ > 0, i?^ > 0. 

For definiteness I choose si > S2. With these parame- 
ters, /3fc-boson is a gapped excitation with gap 

= 6 J(si - S2), (11.13) 

while Qife-boson is the long-range (magnon) excitation 
in the two-spin system. Near the zero wave vector 



E^ 



where the spin-stiffness constant is 

isjJi + 8slJ2 + 2S1S2J 



P = 



Sl 



S2 



(11.14) 



(11.15) 



The spontaneous magnetization of the system M is a sum 
of the spontaneous magnetization on the two sublattices 
M = + , where 



M — < Sf,; > j is from sublattice A 



(11.16) 



< S'lj > j is from sublattic 



B 



In terms of the Holstein-Primakoff bosons (a^ and hk) 
the magnetization adopts the form 

1 



Sl 



N ^ 

k£B, 



-S2 



-E 

keBr 



(11.17) 



< b+bk > 



where Na = Nb = N are the numbers of sites on sub- 
lattices A and B. Finally one can rewrite AI^ and 
in terms of the at (a^ ) and /3fc {P^) excitations 



Sl 



-E 

keBr 



(11.18) 



-S2 



N ^ 

keB, 



^ [ul < 13+ Pk > +vl < aka+ >] 



Then, the magnetization of the system adopts the form 



M ^ si~ S2-j^J2 [< "^"fc > - < ^t^f^ >] 

(11.19) 



keB,^ 

In equations (jII.18III.19P 



< aZak >= 



1 



e^S/T - 1 
1 



< f3tPk >= ^F7^ 



, < aka^ >= 1+ < a^ak > 
(11.20) 

, < PkPt >= 1+ < P+f3k > 



We consider a theory with Hamiltonian Eq. (jl.ip . 
where the sums are over all sites of a three-dimensional 
cubic lattice with space size a = 1. The two equivalent 
sublattices A and B are face centered cubic (fee) lattices 
with space size 2a = 2. The coordinates kx,ky,kz of 
the wave vector k are coordinates in the basis x, y, z, 
which are primitive vectors of the cubic lattice. To im- 
plement the integrations over the wave vector it is more 
convenient to switch to coordinates 91, 92, 93 in the ba- 
sis A, B, C, which arc primitive vectors of the reciprocal 
lattice to the fee lattice with space size equal to 2. 



A — TT {x + y — z) 
B = IT {~x + y + z) 
C — TT {x — y + z) 

To that end we utilize the transformation 

kx = TTqi - TTq2 + 7rgf3 
ky = nqi + TTq2 - nq^ 
kz = -Trgi + 7i'g2 + TT^s 



(11.21) 



(11.22) 



The Jaeobian of the transformation is equal to 47r^, so 
that (fik/{2TT)^ = (fiq/2. When the wave vector k runs 
the reduced Brillouin zone k G Br, the new vector runs 
the cube < qi < 1- Finally we introduce the more 
convenient coordinates pi = irqi; p2 = irq- P3 = 
nq^. They run the interval [0,7r] and (fik/{2TT)^ = 
<Pp/2'K^ . In terms of the new coordinates the fimctions 
of the wave vector (|TL6| and pi.7p adopts the form 

Ep ~ 6 — cos(2pi) — cos(2p2) — cos(2p3) 

- cos2(pi -P2) - cos2(pi -pa) - cos2(p2 - Pz) 

7p = 2 J ^si 32 [cos(-pi +P2+ Pa) + cos(pi - p2 + Pz) 

+ cos(pi +P2 -Ps)] (11.23) 

The magnon excitation-afc in the effective theory 
Eq. (jll.lOp is a complicate mixture of the transversal fluc- 
tuations of the A and B spins. As a result the magnons' 
fluctuations suppress in a different way the magnetic 
order on sublattices A and B. Quantitatively this de- 
pends on the coefficients Uk and Vk in Eqs. (jII.18p . The 
magnetization depends on the dimensionless tempera- 
ture T / J and dimensionless parameters si, S2, Ji/J and 



4 



J2/J. For parameters si = 1.5, S2 = 1, Ji/J = 0.47 
and J2/J = 0.005 the functions M{T/J), M^{T/J) and 
M^{T/J) are depicted in Fig.l. The upper (green) Hne 
is the sublattice A magnetization, the bottom (red) hne 
is the sublattice B magnetization and the middle (blue) 
line is the magnetization of the system. The figure shows 
that the magnetic order on sublattice B (red line) is 
suppressed first, at temperature T* / J = 20. Once sup- 
pressed, the magnetic order can not be restored at tem- 
peratures above T* because of the increasing effect of 
magnon fluctuations. Hence, the sublattice B magneti- 
zation should be zero above T*. It is evident from Fig.l, 
that this is not the result within customary spin-wave 
theory. 



M 



form 




T/ J 



FIG. 1: (color online) Temperature dependence of the ferro- 
magnetic moments: M (blue line)-the magnetization of the 
system, (green line)-sublattice A magnetization, 

(red line)-sublattice B magnetization for parameters si — 
1.5, S2 = 1, Ji/J = 0.47 and J2/J = 0.005: spin-wave 
theory 



III. MODIFIED SPIN- WAVE THEORY 

To solve the problem we use the idea on description 
of paramagnetic phase of 2D ferromagnets (T > 0) by 
means of modified spin- wave theory [3, [l3| • In ths sim- 
plest version the spin-wave theory is modified by intro- 
ducing a new parameter which enforce the magnetization 
of the system to be equal to zero in paramagnetic phase. 

We consider two-sublattice system and to enforce the 
magnetic moments on the two sublattices to be equal to 
zero in paramagnetic pase we introduce two parameters 
Ai and A2. The new Hamiltonian is obtained from the old 
one adding two new terms to the Hamiltonian Eq. ljl.ip 



H 



H - ^Ai5fj + A25'2 



(IH.l) 



-k^-k 



Ik {hak 



keB-r 



where the new dispersions are 
- el + Ai, 



-k 



el 



Utilizing the same transformation Eq 
eters 



A2. 



(in.2) 
(in.3) 

with param- 



Uk 



\ 




(in.4) 



one obtains the Hamiltonian in diagonal forma 



keB,^ 



- 1 



(in.5) 



where 



El 



^{il + 4)2 - A^l - ei + 



{el + eiY ~ 472 + 
{el + elf - 472 



zl 111.6) 



We have to do some assumptions for the parameters Ai 
and A2 to ensure correct definition of the Bose theory. For 
that purpose it is convenient to represent the parameters 
Ai and A2 in the form 

Ai ~ 6Js2/^i — 6Js2, A2 = 6Jsi^2 — 6Jsi. (HI. 7) 



In terms of the new parameters fii and fi2 the dispersions 
; the form 

el = 4siJiefc -I- QJ S2H1 



el and e\ adopt the form 



(in.8 



e,. 



4:S2J2£k + GJ Si^2 



We assume ^1 and ^2 to be positive {fj,i > 0, ^2 > 0), 
then £^ > and e^. > for all values of the wave-vector 
k. The Bose theory is well defined if square-roots in 
equations pil.6|) are well defined and > 0, E^ > 0. 
This comes true if 



^l/i2 > 1. 



(in.9) 



The /3fc-excitation is gapped (i?^ > 0) for all values of 
parameters /ii and fi2 which satisfy Eq. (jIII.9p . The a- 
excitation is gapped if iJ-ifi2 > 1, but in the particular 
case 



In momentum space the new Hamiltonian adopts the 



1 



(IH.IO) 
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Eq = 0; and near the zero wave vector 

« pP (III.ll) 

with spin-stiffness constant 

^ ^ H4J2f^l + gl./lM2) + 2siS2J ^jjj 

In the particular case Eq. ljlll.lOp afc-boson is the long- 
range excitation (magnon) in the system. 

A. si > S2 and Ji » J2 

We introduced the parameters Ai and A2 {ni , /J.2) to en- 
force the spontaneous magnetization on the sublattices A 
and B to be equal to zero. We find out the parameters 
and fi2, as functions of temperature, solving the system 
of two equations = 0, = 

si ~ ^ [^k < alau > +vl < PkPt >] = 

(III.13) 

+ ^ ^ [ul < PkPk > +vl < akal >] = 

where the coefficients and Vk are given by Eqs(|TTL4| 
and the Bose functions pi.20[) are defined with disper- 
sions i;^ and E^. The solutions /ii(T/J) and ^i2{T / J) 
depend on the parameters si, S2, Ji/J and J2I J ■ For 
si = 1.5, S2 = 1, Ji/J = 0.47 and J2/J = 0.005 they are 
depicted in Fig. 2. 




FIG. 2: (color online) The solution of the system of equations 
(|III.13|I iii{T/J), fi2{T/J) and pi{T / J).ijli{T / J) for param- 
eters si = 1.5, S2 = 1, Ji/J = 0.47 and J2/J = 0.005 

The upper (red) line is fi2{T/J), the bottom (black) 
line is fii(T/J) and the middle (blue)line is the product 



fii(T/ J).fi2{T/ J). The numerical calculations show that 
for high enough temperature ^2 > 1, 1 > fJ-i > and 
/ii.jU2 > 1- Hence ak and I3k excitations are gapped. 
When the temperature decreases ^2 decreases remaining 
larger then one, ^1 decreases too remaining positive, and 
the product /ii./i2 decreases remaining larger than one. 
At temperature T/v/J = 50.22 one obtains fii ~ 0.422, 
fi2 = 2.368 and therefor ^1 • ^1 — 1. Hence, at T^r 
long-range excitation (magnon) emerges in the spectrum 
which means that T/v is the Neel temperature. 

Below the Neel temperature the spectrum contains 
magnon excitations, thereupon /ii • /i2 = 1. It is con- 
venient to set 

m = p., ^2 = I/m- (IH.14) 
In ordered phase magnon excitations are origin of sup- 
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FIG. 3: (color online) Temperature dependence of p for T 
between T* and Tat and parameters Si = 1.5, S2 = 1, Ji/J = 
0.47 and J2/J = 0.005 

pression of the magnetization. Near the zero temper- 
ature their contribution is small and at zero tempera- 
ture they are close to si and S2- Increasing the temper- 
ature magnon fluctuations suppress the magnetization. 
For the chosen parameters they first suppress the sublat- 
tice B magnetization at T* {M^{T*) > 0). Once sup- 
pressed, the magnetic moment of sublattice B spins can 
not be restored increasing the temperature above T* . To 
formulate this mathematically, we modify the spin-wave 
theory introducing the parameter /i Eq. (|III.15p . Below 
T* fj, = 1, or in terms of A parameters Ai = A2 = 0, 
which reproduces the customary spin-wave theory. In- 
creasing the temperature above T* the magnetic moment 
of the sublattice B spins should be zero. This is why we 
impose the condition M^{T) = if T > T*. For tem- 
peratures above T* the parameter /x is a solution of this 
equation. The function p,{T / J) is depicted in Fig. 3. 
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Increasing the temperature above T* fj.{T/J) decreases 
from n{T*/J) = 1 to fi{TN/J) = MTn/J) = 0.422. 

Next, one utilizes the so obtained function fi{T/J) to 
calculate the sublattice A magnetization as a function of 
the temperature. Above T* is equal to the mag- 
netization of the system. The magnetic moments of the 
sublattice A and B spins, as well as the magnetization 
of the system, as a function of the temperature arc de- 
picted in Fig. 4 for parameters si = 1.5, S2 = 1, Ji/ J = 
0.47, J2/J = 0.005. 



1,6 
1,4 
1,2 
1,0 

0,. 
0,6 



0,0 
-0,2 
-0,4. 
-0,6 



: MA 


s,=1.5 s^=1 






>^ J/J=047 J2/J=0.005 


■ M 




: MB / 

' 1 ' 





T/J 



FIG. 4: (color online) Temperature dependence of the ordered 
moments: M (blue line)-the magnetization of the system, 
(black line)-sublattice A magnetization, (red line)- 
sublattice B magnetization for parameters si = 1.5, S2 ~ 
1, Ji/J = 0.47 and J2/J = 0.005: modified spin- wave 
theory 



obtain a value close to 2si/iB, where Si is the spin of 
the sublattice A spin operators. The figure shows that 
extrapolations give one and just the same result for all 
values of x. One can accept the fact that the Se con- 
centration do not influence over the value of sublattice A 
spin and si — 1.5. 

Below T* the magnetization is a sum of sublattice A 
and B magnetization. Hence, the magnetization at zero 
temperature is equal to 2(si — 82)^3 ■ Therefore, one 
can determine the sublattice B spin S2- The important 
conclusion is that the effective sublattice B spin S2 de- 
creases with increasing Se concentration. The dimension- 
less magnetization (in Bohr magnetons) per lattice site 
a = 2{M^ + M^) is calculated for si = 1.5, Ji/J = 
0.47, J2/J = 0.005 and sa = 1; 0.7; 0.4. The curves 
are depicted in Fig. 5. The figure shows that the present 



a 
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To compare the theoretical results and the experi- 
mental M{T) curves one has, first of all, to interpret 
adequately the measurements. The magnetic moments 
in some materials are close to "spin only" value 2^bS 
and the sublattice spins si and S2 can be obtained 
from the experimental curves. I consider the system 
MnCr2Si-xSex. It has been investigated by measure- 
ments of the magnetization at 15.3A:Oe as a function of 
temperature (figure 94 in [^). The maximum in the mag- 
netization versus temperature curve, which is typical of 
MnCr2S4 (x = 0), increases when x increase, and disap- 
pears at X — 0.5. The Neel temperature decreases from 
74i^ at a; = to 56K at x = 2. The authors' conclusion 
is that the observed change of the magnetic properties is 
attributed to a decrease of the strength of the negative 
Mn^"*" — Cr"^"*" superexchange interaction with increasing 
Se concentration. 

As follows from the present theory (see Fig. 4 middle- 
blue line) the maximum of the magnetization is at T*. 
Above T* the magnetization of the system is equal to 
the magnetization of sublattice A spins. If we extrapo- 
late this curve below T* up to zero temperature we will 



FIG. 5: (color online)The magnetization, in Bohr magne- 
tons, per lattice site a = 2{M'^ + M^) for si = 1.5, Ji/J = 
0.47, J2/J = 0.005 and sa = 1; 0.7; 0.4 

calculations capture the essential features of the system; 
increasing 5*6 concentration (decreasing S2) leads to de- 
creasing of Neel temperature, T* temperature decreases 
too, and the maximum of the magnetization a(T*) in- 
creases. 



B. si > S2 and J2 2> Ji 

Next I consider the case when the sublattice A spin 
si is larger then the sublattice B spin 82 but the intra- 
sublattice exchange constant Ji is much smaller than 
the intra-sublattice B exchange constant J2. Increasing 
the temperature, the magnon fluctuations suppress the 
magnetic order and sublattice A magnetization decreases 
faster than sublattice B magnetization. There is a tem- 
perature Tc at which the magnetization of the system 
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is zero M{Tc) = (the compensation point). Increasing 
the temperature above the sublattice A magnetization 
becomes equal to zero at T*. Above this temperature the 
sublattice A magnetization should be kept equal to zero 
and one utilizes the modified spin-wave theory to calcu- 
late the magnetization of the system which is equal to 
sublattice B magnetization. The magnetization curves 
M^{T/J), M^{T/J) and M{T/J) are depicted in Fig. 
6 for parameters si = 18, S2 = 3, ji = Ji/J = 0.005 and 

J2 = J2/J=l. 




-4-1 • ^ • 1 • ^ • 1 • 1 • ^ • 1 • 1 T/J 

100 200 300 400 500 600 700 800 



FIG. 6: (color online) The magnetization curves: M'* (T/J) - 
upper (black) line, (T/J) - bottom (red) line, M{T/J)- 
middle (blue) line for parameters si — 18, S2 = 3, ji — 
Ji/J = 0.005 and j2 = J2/J = 1: modified spin-wave 
theory 

To compare the theoretical results with experimental 
ones I address the magnetization-temperature curves of 
rare earth iron garnets. When a rare earth ion, such as 
Gd^^ through Yb^^, is present instead of a diamag- 
netic ion, such as y3+ Yttrium Iron Garnet, it is 
found that the magnetization as a function of temper- 
ature shows compensation points, i.e., temperatures at 
which the spontaneous magnetization is zero (Fig. 5 
I use the "spin only" value for the magnetization, and the 
sublattices' spins are consider as an effective parameters 
to be determined from the experimental curves. Usu- 
ally the magnetization, in Bohr magnetons, per lattice 
site a = 2\M^ + M^\ are depicted in figures. T* is 
the temperature at which the sublattice A magnetiza- 
tion becomes equal to zero (black line Fig. 6), and hence 
the magnetization is minimal (blue line Fig. 6), in the 
case of the cr(T') curves T* is the temperature, above the 
compensation point, at which a(T*) is maximal. The 
magnetization above T* is equal to the sublattice B mag- 
netization. To assess the value of the sublattice B spin 
S2, one has to extrapolate the curve to zero tempera- 
ture. Then the spin si can be obtained from the zero 
temperature of the magnetization 2(si — 82)^.3- The ex- 



perimental curves (Fig. 5 [4]) show that the value of 
a{T*) is different for different rare earth ions, and hence 
the effective spin §2 is different for different rare earths. 
It is smaller for gadolinium, increases for terbium, dys- 
prosium, holmium, and is higher for ytterbium. At the 
same time the temperature T* decreases, which means 
that sublattice B intra-cxchange constant J2 decreases. 
The magnetization a — 2\M^ + for two choices 

of the parameters are depicted in Fig. 7: black squares- 
si = 18, S2 = 3, ji = Ji/J = 0.005, j2 = J2/J = 1 and 
red circles-si = 18, S2 = 6, ji = Ji/J = —0.05, j2 = 
J2/J = 0.15. The negative sign of Ji in the second 
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FIG. 7: (color online)The magnetization a = 2\M^ + M'^\ as 
a function of T/J: black squares- si = 18, S2 = 3, Ji/J = 
0.005, J2/J = 1, red circles -si = 18, S2 = 6, Ji/J = 
-0.05, J2/J = 0.15 

group of parameters is chosen to reproduce the shape 
of the experimental curves. One can do better fitting, 
with positive constants J2 for all rare earths, accounting 
for magnon scattering processes. 



C. si — S2 and Ji 3> J2 

Finally I consider a theory with parameters si = S2 = 
s and Ji S> J2- It is more correctly to thought of this 
case as a generalization of the antiferromagnetism. The 
system has two magnons with dispersions and E^. 
Near the zero wave vector the dispersions' asymptotic is 
as in the antiferromagnetic case 

E^^vs\kl E^^^vs\k\ (III.15) 
with spin-wave velocity 



Vs = 2sv/l2Ji + I2J2 + 3J. (III.16) 
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As a result, at zero temperature and near the zero tem- 
perature one observes a compensation of sublattice A and 
B magnetization and magnetization of the system is zero. 
The difference of the dispersions is 



4s(Ji - J2)£fc « (Ji - J2)k2, (III.17) 



therefore increasing the temperature one obtains non- 
compensation of sublattice A and B magnetization. The 
magnon fluctuations suppress first the sublattice B mag- 
netization at temperature T*. Above this temperature 
only sublattice A contributes the magnetization. The 
magnetization curves M^{T/J), M^{T/J) and M{T/ J) 
are depicted in Fig. 8 for parameters si = 1, S2 = 
1, Ji/J = 0.47 and J2/J = 0.005. The magnetization- 
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FIG. 8: (color online) The magnetization curves: M'*(r/J) - 
black squares, (T/J) -red circles, M(r/J)-blue triangles 
for parameters si = 1, S2 = 1, Ji/J = 0.47 and J2/J = 0.005: 
modified spin- wave tiieory 



temperature curve (blue triangles) reproduces the Neel's 
anomaly. 



IV. SUMMARY 

In summary, I have calculated the magnetization as a 
function of temperature for two-sublattice ferrimagnet. 
Anomalous temperature dependence of the magnetiza- 
tion, predicted by Ncel, is reproduced, and the shape 
of the theoretical curves satisfactorily coincides with the 
experimental ones. The most important difference be- 
tween Neel's theory and the present modified spin-wave 
theory is that Neel's calculations predict a temperature 
T/v at which both the sublattice A and B magnetiza- 
tions become equal to zero. The modified spin wave 
theory predicts to phases: at low temperatures (0,T*) 
the magnetic orders of the two sublattices contribute to 



the magnetization of the system, while at high temper- 
atures (r*,T/v) only one of the sublattices, that with 
stronger intra-sublattice exchange, has non-zero sponta- 
neous magnetization. It is important to stress that de- 
spite the fact that only one of the spins contribute the 
magnetic order at high temperatures T* < T < Tn this is 
not a ferromagnetic phase. This is a ferrimagnctic phase 
because the magnon of the system is a mixture of the 
transversal fiuctuations of the two spins. 

Two ferromagnetic phases where theoretically pre- 
dicted, very recently, in spin-Fermion systems, which 
obtain their magnetic properties from a system of lo- 
calized magnetic moments being coupled to conducting 



electrons . At the characteristic temperature T* , the 
magnetization of itinerant electrons becomes zero, and 
high temperature ferromagnetic phase (T* <T <Tc) is 
a phase where only localized electrons give contribution 
to the magnetization of the system. An anomalous in- 
creasing of magnetization below T* is obtained in good 
agrement with experimental measurements of the ferro- 
magnetic phase of [/G'e2[l2|- The results of the present 
paper and the previous one [ll| suggest that T* transi- 
tion from a magnetic phase to another magnetic phase is 
a generic feature of the two spin systems. The additional 
phase transition demonstrates itself through the anoma- 
lous temperature variation of the spontaneous magneti- 
zation. Another way is possible, to experimentally iden- 
tify the transition, if T* is a temperature at which the 
itinerant electrons in the system start to form magnetic 
moment. We know from itinerant ferromagnetism that 
the transition is observed through the change in the tem- 
perature dependence of resistivity. The most prominent 
example is the spinel Fe^Oi. It has been extensively in- 
vestigated and most striking feature is the Verwcy transi- 
tion 111]. At relatively low temperature Tt = 100- 
the magnetization abruptly changes the slope, and the 
conductivity has anomalous behavior. Many decades of 
research on FesO^ have led to the view [IJ] that the 
spinel is a two-sublattice system with three spins. The 
Fe^^ and Fe^^ ions are associated with localized spins 
on sublattices A and B, while Fe^^^ ions are associated 
with itinerant electrons on sublattice B. With this in 
mind one can interpret Tf as a temperature at which the 
magnetic moment of itinerant electrons sets in. 

The present theory of ferrimagnetism permits to for- 
mulate phenomenological rules which to enable the anal- 
ysis of experimental results for more complicate systems 
without explicit calculations. For example, the three- 
spin systems have, most generally, three characteristic 
temperatures T* < T2 < Tn- The Neel temperature Tn 
is the temperature at which the magnetic order sets in. 
When T2 < T < Tjy only one of the spins has non zero 
spontaneous magnetization and the magnetization of the 
system equals the magnetic order of that spin. When 
T* < T < T2 two of the spins contribute the magne- 
tization of the system. If AI{T2) is local maximum the 
two spins are with anti-parallel order. If the magneti- 
zation shows an anomalous enhancement below T2 the 
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two spins are parallel. Below Tj* all three spins have non 
zero spontaneous magnetization, and the magnetization 
changes the slope at this temperature. 

As an example I consider the CeCrSb^ compound [lH . 
The magnetic behavior is determined by Cr ions and 
Ce — Af electrons, which are sat in different sublattices. 
The measurements of the magnetization and electric re- 
sistivity on single crystal of CeCrSb^ show that there 
arc three characteristic temperatures, Neel temperature 
Tn = UbK, T* = 108K and T* = 18K. The close- 
ness of T/v and T2 leads to misinterpretation and even to 
identification of these temperatures, but evident increase 
of magnetization when T2 < T < Tn and subsequent 
decrease below Tj* proves that there are three character- 
istic temperatures. Below T* the magnetization increases 
again. The temperature variation of the resistivity shows 
anomalous behavior at and rj* . The experimental re- 
sults suggest that part of the itinerant Cr electrons start 
to form magnetic order at Neel temperature Tjv and the 
rest one at Tj*, while the locahzed Ce — 4/ electrons do 
this at T2*. At high temperature (T2*,Tjv) the magnetic 
order of some of the Cr electrons contribute the mag- 
netization of the system. Below Tj* the the magnetic 
moment of localized Ce electrons sets in. The low tem- 



perature phase (T*,T2*) is a phase where all Cr and Ce 
electrons contribute the magnetization of the system. 

Another example is the manganese vanadium oxide 
spinel MnV20i ^M,- C)n one of the sublattices sites is 
Afn+^ in a Sd^ high spin configuration s = 5/2 with 
quenched orbital angular momentum. On the sites of 
the other sublattices are the vanadium electrons. More 
precise measurements of the magnetization as a function 
of temperature [13] show that the set in of magnetic or- 
der is due to Mn localized electrons at Neel temperature 
Tjv = b&.bK. Because of the splitting of the vanadium 
three bands some of electrons start to form magnetic mo- 
ment at T2* = A%K . An evidence for this is the abrupt 
decrease of magnetization below Tj*, which also indicates 
that the magnetic order of vanadium electrons is anti- 
parallel with the order of Mn electrons. Below T* = 27 K 
the magnetization decreases faster which means that all 
vanadium electrons are involved in the magnetization of 
the system. 
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